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Abstract 


Nomenclature 


Experimental  data  taken  from  gas  turbine  engines  has 
shown  that  hot  streaks  exiting  combustors  can  have  a 
significant  impact  upon  the  secondary  flow  and  wall 
temperature  of  the  first  stage  turbine  rotor.  Under¬ 
standing  the  secondary  flow  and  heat  transfer  effects 
due  to  combustor  hot  streaks  is  essential  to  turbine 
designers  attempting  to  optimize  turbine  cooling  sys¬ 
tems.  A  numerical  investigation  is  presented  which  ad¬ 
dresses  the  issues  of  multi-blade  count  ratio  and  three- 
dimensionality  effects  on  the  prediction  of  combustor 
hot  streak  migration  in  a  turbine  stagey  The  two-  and 
three-dimensional  Navier-Stokes  analyses  developed  by 
Rai  et  al  are  used  to  predict  unsteady  viscous  rotor- 
stator  interacting  flow  in  the  presence  of  a  combus¬ 
tor  hot  streak.  Predicted  results  are  presented  for 
a  two-dimensional  3-stator/4-rotor,  a  two-dimensional 
1-stator/l-rotor,  and  a  three-dimensional  1-stator/l- 
rotor  simulation  of  hot  streak  migration  through  a  tur¬ 
bine  stage.  Comparison  of  these  results  with  experi¬ 
mental  data  demonstrates  the  capability  of  the  three- 
dimensional  procedure  to  capture  most  of  the  flow 
physics  associated  with  hot  streak  migration  including 
the  effects  of  combustor  hot  streaks  on  turbine  rotor 
surface  temperatures. 
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-  Speed  of  sound 

-  Specific  energy 

-  Total  energy 

-  Mach  number 

-  Static  pressure 

-  Prandtl  number 

-  Stagnation  pressure 

-  Free  stream  inlet  reference  Reynolds  number 

-  Static  temperature 

-  x  component  of  velocity 

-  Rotor  velocity 

-  y  component  of  velocity 

-  z  component  of  velocity 

-  thermal  conductivity 

-  Second  coefficient  of  viscosity 

-  First  coefficient  of  viscosity 

-  Density 

-  Shear  stress 

-  Rotor  rotational  speed 
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-  Hot  streak 

-  Inviscid 

-  Laminar  quantity 

-  Stagnation  quantity 

-  Turbulent  quantity 

-  Viscous 

-  First  derivative  with  respect  to  x,  y  or  z 

-  Second  derivative  with  respect  to  x,  y  or  z 

-  Inlet  quantity 

-  Exit  quantity 
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Introduction 

The  drive  toward  low  fuel  consumption  and  maxi¬ 
mum  gas  turbine  efficiency  for  commercial  aircraft  has 
led  to  increases  in  combustor  exit  temperature  which  in 
turn  has  had  a  direct  impact  on  the  durability  of  first 
stage  turbine  airfoils.  In  order  to  keep  turbine  airfoil 
wall  temperatures  below  critical  levels,  complex  inter¬ 
nal  cooling  schemes  are  designed  to  lay  a  film  of  rela¬ 
tively  cool  air  adjacent  to  the  blade.  This  film  cooling 
air  acts  as  a  sealant  to  keep  the  high  temperature  com¬ 
bustor  flow  from  coming  in  direct  contact  with  turbine 
airfoil  surfaces.  During  the  design  of  the  first  stage  tur¬ 
bine,  flow  field  analyses  which  assume  a  uniform  turbine 
inlet  temperature  profile  are  used  to  predict  airfoil  sur¬ 
face  temperature  distributions.  These  wall  temperature 
distributions  are  subsequently  used  to  guide  the  design 
of  internal  cooling  schemes  in  terms  of  defining  cooling 
mass  flow  requirements  and  its  distribution. 

Experimental  data  taken  from  actual  gas  turbine 
combustors,  however,  indicate  that  the  flow  exiting  the 
combustor  has  both  circumferential  and  radial  tempera¬ 
ture  gradients.  These  temperature  gradients  arise  from 
the  combination  of  the  combustor  core  flow  with  the 
combustor  bypass  and  combustor  surface  cooling  flows. 
It  has  been  shown  both  experimentally  and  numeri¬ 
cally  [1,2]  that  temperature  gradients,  in  the  absence 
of  total  pressure  non-uniformities,  do  not  alter  the  flow 
within  the  first  stage  turbine  stator  but  do  have  signifi¬ 
cant  impact  on  the  secondary  flow  and  wall  temperature 
of  the  first  stage  rotor.  Combustor  hot  streaks,  which 
can  typically  have  temperatures  twice  the  free  stream 
stagnation  temperature,  increase  the  extent  of  the  sec¬ 
ondary  flow  in  the  first  stage  rotor  and  significantly  alter 
the  rotor  surface  temperature  distribution.  A  combus¬ 
tor  hot  streak  such  as  this  has  a  greater  streamwise  ve¬ 
locity  than  the  surrounding  fluid  and  therefore  a  larger 
positive  incidence  angle  to  the  rotor  as  compared  to 
the  free  stream.  Due  to  this  rotor  incidence  variation 
through  the  hot  streak  and  the  slow  convection  speed 
on  the  pressure  side  of  the  rotor,  the  hot  streak  accu¬ 
mulates  on  the  rotor  pressure  surface  creating  a  “hot 


can  be  useful  for  both  increasing  the  understanding  of 
this  phenomena  as  well  as  to  guide  engineers  for  fu¬ 
ture  experimental  investigations.  The  focus  of  the  cur¬ 
rent  effort  has  been  to  perform  both  two-  and  three- 
dimensional  rotor/stator  interaction  numerical  simula¬ 
tions  which  include  the  effects  of  combustor  hot  streaks 
in  order  to  provide  additional  insight  into  the  hot  streak 
migration  phenomena. 

In  the  present  investigation,  the  two-  and  three- 
dimensional  unsteady  Navier-Stokes  analyses  of  Rai 
et  al  [3,4]  are  used  to  predict  the  effects  of  a  com¬ 
bustor  hot  streak  on  the  flow  in  the  first  stage  of 
a  turbine.  Two-dimensional  simulations,  including  3- 
stator/4-rotor  and  1-stator/l-rotor  configurations,  have 
been  performed  to  study  the  effects  of  stator/rotor  blade 
count  on  the  time-averaged  rotor  surface  mid-span  tem¬ 
perature  distribution.  In  addition,  to  determine  the 
three-dimensional  effects  of  the  hot  streak  migration 
on  the  secondary  flow  and  rotor  surface  temperature 
distribution,  a  three-dimensioned  i-stator/l-rotor/l-hot 
streak  simulation  has  been  performed.  The  results  from 
these  computations  have  been  compared  to  the  experi¬ 
mental  data  reported  by  Butler  et  al  [1]  and  the  compu¬ 
tations  performed  by  Rai  and  Dring  [2],  Krouthen  and 
Giles  [5],  and  Takahashi  and  Ni  [6,7].  These  numerical 
results  and  their  comparison  with  the  available  exper¬ 
imental  data  sheds  new  insight  into  the  phenomena  of 
hot  streak  migration  and  its  effect  on  first  stage  turbine 
aerodynamics  and  heat  transfer. 


Numerical  Integration  Procedure 

The  governing  equations  considered  in  this  study  are 
the  time  dependent,  three-dimensional  Navier-Stokes 
equations: 


U,  +  (Fi  +  Fv)x  +  (Gi  +  Gv)v  +  (Hi  +  Hv)t=  0  (1) 


spot”. 

Both  the  secondary  flow  and  wall  temperature  shift-  where 
ing  effects  due  to  combustor  hot  streaks  are  important 
phenomena  which  should  be  accounted  for  by  turbine 
designers  to  increase  turbine  aerodynamic  performance 
and  optimize  turbine  cooling  flow  schemes.  Without 
properly  accounting  for  non-uniformities  in  turbine  inlet 
temperature  profiles  such  as  those  related  to  combustor 
hot  streaks,  first  stage  turbine  performance  and  dura¬ 
bility  could  suffer.  Computational  fluid  dynamic  sim¬ 
ulations  of  hot  streak  migration  in  turbine  stage  flows 
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For  the  present  application,  the  second  coefficient 
of  viscosity  is  calculated  using  Stokes’  hypothesis, 
A  =  -2/3p.  The  equations  of  motion  are  completed  by 
the  perfect  gas  law. 

The  viscous  fluxes  are  simplified  by  incorporating  the 
thin  layer  assumption  [8].  In  the  current  study,  viscous 


terms  are  retained  in  the  direction  normal  to  the  hub 
surface  (z  direction)  and  in  the  direction  normal  to  the 
blade  surface  (y  direction).  To  extend  the  equations  of 
motion  to  turbulent  flows,  an  eddy  viscosity  formula¬ 
tion  is  used.  Thus,  the  effective  viscosity  and  effective 
thermal  conductivity  can  be  defined  as: 


n  =  HL+  HT 


(7) 


The  turbulent  viscosity,  /ij,  is  calculated  using  the 
Baldwin-Lomax  [8]  algebraic  turbulence  model  as  ap¬ 
plied  by  Hung  and  Buning  [9]  for  three-dimensional 
flows. 

The  numerical  procedure  for  the  three-dimensional 
analysis  consists  of  a  time  marching,  implicit,  third- 
order  spatially  accurate,  upwind,  finite  difference 
scheme.  The  inviscid  fluxes  are  discretized  using  a  com¬ 
bination  of  central,  forward,  and  backward  differences 
depending  on  the  local  eigenvalues  of  the  flux  Jaco- 
bians  according  to  the  scheme  developed  by  Roe  [10]. 
The  viscous  fluxes  are  calculated  using  standard  cen¬ 
tral  differences.  The  alternate  direction,  approximate- 
factorization  technique  of  Beam  and  Warming  [11]  is 
used  to  compute  the  time  rate  changes  in  the  primary 
variables.  An  inner  Newton  iteration  can  be  used  to  in¬ 
crease  stability  and  reduce  linearization  errors.  For  all 
cases  investigated  in  this  study,  three  Newton  iterations 
were  performed  at  each  time  step.  The  numerical  pro¬ 
cedure  used  for  the  two-dimensional  analysis  is  similar 
to  that  used  in  the  three-dimensional  analysis,  except 
that  the  inviscid  fluxes  are  discretized  according  to  the 
scheme  developed  by  Osher  [12].  Further  details  of  the 
two-  and  three-dimensional  numerical  techniques  can  be 
found  in  Refs.  [2,3,4]. 


Grid  Generation  and  Geometry 

The  three-dimensional  Navier-Stokes  analysis  devel¬ 
oped  by  Rai  et  al  [2,3,4]  uses  five  zonal  grids  to  discretize 
the  rotor-stator  flow  field  and  facilitate  relative  motion 
of  the  rotor.  A  combination  of  O-  and  H-grid  sections 
are  generated  at  constant  radial  spanwise  locations  in 
the  blade-to-blade  direction  extending  upstream  of  the 
stator  leading  edge  to  downstream  of  the  rotor  trailing 
edge.  Algebraically  generated  H-grids  are  used  in  the 
regions  upstream  of  the  leading  edge,  downstream  of 
the  trailing  edge  and  in  the  inter-blade  region.  The  O- 
grids,  which  are  body-fitted  to  the  surfaces  of  the  air¬ 
foils  and  generated  using  an  elliptic  equation  solution 
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procedure,  are  used  to  properly  resolve  the  viscous  flow 
in  the  blade  passages  and  to  easily  apply  the  algebraic 
turbulence  model.  Computational  grid  lines  within  the 
O-grids  are  stretched  in  the  blade-normal  direction  with 
a  fine  grid  spacing  at  the  wall.  The  combined  H-  and 
O-  overlaid  grid  sections  are  stretched  in  the  spanwise 
direction  away  from  the  hub  and  tip  regions  with  a  fine 
grid  spacing  located  adjacent  to  the  hub  and  tip.  An 
O-grid  structure  is  used  in  the  region  between  the  ro¬ 
tor  blade  and  the  tip  end  wall  in  the  tip  clearance  region. 
Grid  generation  for  the  two-dimensional  analysis  is  sim¬ 
ilar  to  that  for  the  three-dimensional  procedure  except 
that  the  O-grids  are  patched,  rather  than  overlaid,  into 
the  H-grids. 

Boundary  Conditions 

The  theory  of  characteristics  is  used  to  determine  the 
boundary  conditions  at  the  stator  inlet  and  rotor  exit. 
For  subsonic  inlet  flow,  the  total  pressure,  v  and  w 
velocity  components  and  the  downstream  running  Rie- 
mann  invariant,  Ri  =  u  +  are  specified  while  the 
upstream  running  Riemann  invariant,  R?  =  u  — 
is  extrapolated  from  the  interior  of  the  computational 
domain.  Inlet  flow  boundary  conditions  within  the 
hot  streak  region  are  updated  using  Riemann  invari¬ 
ants  where  the  velocity  and  speed  of  sound  reflect  the 
static  and  stagnation  temperature  increase  in  the  hot 
streak.  The  static  and  total  pressure  in  the  hot  streak 
is  assumed  to  be  equal  to  that  of  the  undisturbed  inlet 
flow.  For  subsonic  outflow  the  pressure  ratio  Pi/Pti  is 
specified  while  the  v  and  w  velocity  components,  en¬ 
tropy,  and  the  downstream  running  Riemann  invariant 
are  extrapolated  from  the  interior  of  the  computational 
domain.  Absolute  no-s!ip  boundary  conditions  are  en¬ 
forced  at  the  hub  and  tip  endwalls  of  the  stator  region, 
along  the  surface  of  the  stator,  and  along  the  outer  cas¬ 
ing  (tip  endwall)  of  the  rotor.  Relative  no-slip  bound¬ 
ary  conditions  are  imposed  at  the  hub  and  along  the 
surface  of  the  rotor.  Periodicity  is  enforced  along  the 
outer  boundaries  of  the  H-grids  in  the  circumferential 
(0)  direction.  In  the  present,  study  the  flow  is  assumed 
to  be  adiabatic. 

Dirichlet  conditions,  in  which  the  time  rate  change  in 
the  vector  U  of  Eq.  (2)  is  set  to  zero,  are  imposed  at 
the  grid  interface  boundaries  located  between  the  stator 
and  rotor  H-grids  and  at  the  overlaid  boundaries  of  the 
O-  and  H-grids  in  the  stator  and  rotor  regions.  The  flow 
variables  of  U  at  zonal  boundaries  are  explicity  updated 
after  each  time  step  by  interpolating  values  from  the  ad¬ 
jacent  grid.  Because  of  the  explicit  application  of  the 
zonal  boundary  conditions,  large  time  steps  necessitate 


the  use  of  more  than  one  Newton  iteration.  The  zonal 
boundary  conditions  are  non-conservative,  but  for  sub¬ 
sonic  flow  this  should  not  affect  the  accuracy  of  the  final 
flow  solution.  Further  information  describing  the  imple¬ 
mentation  of  the  boundary  conditions  can  be  found  in 
Refs.  [2,3,4,13]. 

Results 

A  numerical  investigation  of  hot  streak  migration 
has  been  conducted  in  which  predicted  results  using 
the  two-  and  three-dimensional  versions  of  the  current 
Navier-Stokes  procedure  (ROTOR2,  ROTOR3)  have 
been  compared  with  experimental  data  reported  by 
Butler  et  al  [1].  The  turbine  airfoil  geometry  used  in 
both  the  two-  and  three-dimensional  analyses  consists 
of  the  first  stage  of  the  UTRC  Large  Scale  Rotating  Rig 
(LSRR)  [1,14]  which  includes  22  stator  airfoils  and  28 
rotor  airfoils.  An  accurate  simulation  of  this  configura¬ 
tion  requires  at  least  11  stator  and  14  rotor  airfoils.  For 
the  1-stator/l-rotor  two-  and  three-dimensional  simu¬ 
lations,  the  rescaling  strategy  of  Madavan  et  al  [15]  is 
used  to  reduce  the  number  of  airfoils  to  one  stator  and 
one  rotor.  The  stator  was  scaled  down  by  a  factor  of 
(22/28)  and  it  was  assumed  that  there  were  28  stator 
airfoils.  The  pitch  to  chord  ratio  of  the  stator  was  not 
changed.  This  is  different  than  the  scaling  strategy  used 
by  Rai  and  Dring  [2]  in  which  the  rotor  was  scaled  by 
a  factor  of  (28/22)  and  it  was  assumed  there  were  22 
rotor  airfoils.  For  the  two-dimensional  3-stator/4-rotor 
simulation,  a  rescaling  strategy  was  used  to  reduce  the 
number  of  airfoils  to  3  stators  and  4  rotors.  In  this  case 
it  was  assumed  that  there  were  21  stator  airfoils  and  28 
rotor  airfoils  and  the  stators  were  enlarged  by  the  factor 
(22/21). 

In  the  experimental  study  [1],  one  hot  streak  was 
introduced  through  a  circular  pipe  at  40%  span,  and 
midway  between,  two  stator  airfoils  of  the  LSRR.  The 
temperature  of  the  hot  streak  was  twice  that  of  the  sur¬ 
rounding  inlet  flow,  whereas  th>  hot  streak  static  and 
stagnation  pressures  were  identical  to  the  free  stream. 
The  hot  streak  was  seeded  with  CO2  and  the  path  of 
the  hot  streak  determined  by  measuring  CO2  concen¬ 
trations  at  various  locations  within  the  turbine  stage. 
In  the  two-dimensional  1-stator/l-rotor  simulation  and 
three-dimensional  1-stator/l-rotor  simulation,  one  hot 
streak  is  introduced  to  the  inlet  of  each  stator  pas¬ 
sage.  -For  the  two-dimensional  3-stator /4-rotor  simu¬ 
lation,  one  hot  streak  is  introduced  to  the  inlet  of  every 
third  stator  passage.  The  temperature  of  the  hot  streak 
in  the  current  investigation  is  1 .2  times  the  temperature 
of  the  surrounding  inlet  flow.  For  the  two-dimensional 
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3-stator/4-rotor  simulation,  the  hot  streak  was  intro¬ 
duced  over  a  distance  equal  to  one  quarter  of  a  stator 
pitch  and  centered  at  the  mid-gap  position  between  the 
second  and  third  stator.  For  the  two-dimensional  1- 
stator/l-rotor  simulation,  the  hot  streak  was  introduced 
over  one  quarter  of  the  stator  pitch  and  centered  at  mid¬ 
gap.  In  the  three-dimensional  simulation,  the  hot  streak 
was  introduced  through  a  two  inch  diameter  circular  re¬ 
gion  which  corresponds  to  the  two  inch  diameter  pipe 
used  in  the  experiment  [1],  The  center  of  the  hot  streak 
was  located  at  the  mid-gap,  40%  span  location.  In  both 
the  two-  and  three-dimensional  calculations,  the  tem¬ 
perature  profile  used  at  the  stator  inlet  to  simulate  the 
combustor  hot  streak  consisted  of  a  hyperbolic  tangent 
(step-like)  distribution. 

A  15%  axial  gap  between  the  stator  and  rotor  was 
used  in  both  the  two-dimensional  and  three-dimensional 
simulations,  while  the  experimental  configuration  had  a 
65%  axial  gap.  It  has  been  demonstrated  [16]  that  the 
axial  gap  has  no  significant  impact  on  either  the  time- 
averaged  pressure  distributions  or  time-averaged  heat 
transfer  coefficients  on  the  stator  or  rotor.  Secondary 
flow  and  other  viscous  mechanisms,  however,  can  be 
affected  by  the  stator/rotor  axial  gap.  The  inlet  Mach 
number  to  the  stator  was  0.07  and  the  inlet  flow  was 
assumed  to  be  axial.  The  rotor  rotational  speed  was 
410  rpm.  The  experimental  flow  coefficient  was  4>  = 
ti/CI  =  .68,  while  the  numerical  simulations  held  a  flow 
coefficient  of  <j>  =  .78.  The  free  stream  Reynolds  number 
was  100,000  per  inch.  A  pressure  ratio  of  Pi/Pti  = 
.963  was  determined  from  the  inlet  total  pressure  and 
the  static  pressure  measured  in  the  rotor  trailing-edge 
plane. 

Two-Dimensional  Simulations 

A  two-dimensional  3-stator/4-rotor  and  a  1-stator/l- 
rotor  hot  streak  simulation  have  been  performed  as 
part  of  a  numerical  investigation  into  the  effects  of  sta¬ 
tor  to  rotor  blade  count  on  hot  streak  migration  and 
time-averaged  wall  temperature  distributions.  The  two- 
dimensional  calculations  performed  during  this  investi¬ 
gation  were  computed  on  a  four  processor  Aliiant  FX- 
80  mini-supercomputer.  Typical  calculations  required 
.00191  seconds  per  grid  point  per  iteration  computation 
time.  For  both  two-dimensional  simulations,  approxi¬ 
mately  six  cycles  at  3000  iterations  per  global  cycle  were 
needed  to  obtain  a  time-periodic  solution.  A  global  cy¬ 
cle  corresponds  to  the  rotor  blade  rotating  through  "an 
angle  of  2 x/N  where  N  is  the  number  of  stator  blades 
(i.e.  N=3  and  N=l). 

For  the  two-dimensional  3-stator/4-rotor  simulation, 


each  of  the  three  stator  grids  was  constructed  with  101  x 
31  (streamwise  x  tangential)  grid  points  in  the  O-grid 
and  71  x  51  grid  points  in  the  H-grid.  Each  of  the  four 
rotor  grids  was  constructed  with  101  x  31  grid  points  in 
the  O-grid  and  75  x  51  grid  points  in  the  H-grid.  A  total 
of  48,080  grid  points  and  an  average  blade  wall  spacing 
of  2.0  x  10-4  inches  were  used  in  this  simulation.  For  the 
two-dimensioned  1-stator/l-rotor  simulation,  the  stator 
zone  was  constructed  with  101x31  grid  points  in  the  O- 
grid  and  75x31  grid  points  in  the  H-grid.  The  rotor  zone 
was  constructed  with  101  x  31  grid  points  in  the  O-grid 
and  71  x  31  grid  points  in  the  H-grid.  A  total  of  10,788 
grid  points  and  an  average  wall  spacing  of  2.0  x  10~4 
inches  were  used  in  this  calculation.  Examples  of  the 
grid  topology  used  in  the  two-dimensional  simulations 
are  shown  in  Refs.  [2,4]. 

Results  of  the  two-dimensional  simulations  have 
been  compared  with  the  time-averaged  experimental 
data  [1,14].  Figure  1  shows  a  comparison  of  the  time- 
averaged  stator  surface  pressure  distribution  predicted 
by  the  numerical  analysis  for  the  3-stator/4- rotor  and 
1-stator/l-rotor  configurations  with  experimental  data 
reported  by  Dring  et  al  [14].  The  time-averaged  pres¬ 
sure  coefficient  is  defined  as: 


Cp  — 


P avg  ~  Pfl 

hPxU* 


(8) 


where  Pavg  if  the  local  time-averaged  pressure,  Pt\  is 
the  inlet  total  pressure,  p\  is  the  average  inlet  free 
stream  density,  and  U  is  the  rotor  velocity.  The  pre¬ 
dicted  results  for  the  3-stator /4-rotor  and  1-stator/l- 
rotor  configurations  both  exhibit  good  agreement  with 
the  experimental  data.  Figure  2  presents  a  compari¬ 
son  of  the  predicted  results  of  the  two  simulations  with 
the  experimental  time-averaged  surface  pressure  data 
for  the  rotor.  Excellent  agreement  exists  between  the 
numerical  predictions  and  the  data.  The  predicted  re¬ 
sults  shown  in  Figs.  1  and  2  are  also  in  good  agreement 
with  those  reported  by  Rai  et  al  [2,4]  for  a  similar 
1-stator/l-rotor  configuration  with  and  without  a  hot 
streak.  These  figures  substantiate  that  the  effect  of  the 
hot  streak  on  the  time-averaged  pressure  fields  of  the 
stator  and  rotor  is  negligible  since  the  static  and  stag¬ 
nation  pressure  of  the  hot  streak  is  identical  to  that  of 
the  free  stream. 

A  measure  of  the  unsteadiness  of  the  flow  can  be  ob¬ 
tained  by  evaluating  the  size  of  the  surface  pressure  flu- 
cuations  of  the  stator  and  rotor.  The  pressure  fluctua¬ 
tions  can  be  quantified  through  the  use  of  an  unsteady 
pressure  amplitude  coefficient  defined  by: 
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Figure  1:  Predicted  and  experimental  blade  loading  for 
the  stator  (2-D) 


Figure  2:  Predicted  and  experimental  blade  loading  for 
the  rotor  (2-D) 


Figure  3:  Pressure  amplitude  coefficient  distribution  for 
stator  (2-D) 


Figure  4:  Pressure  amplitude  coefficient  distribution  for 
rotor  (2-D) 


experimental  data.  This  figure  shows  that  the  ampli- 
p  -p  tude  of  the  pressure  variations  is  small  on  the  stator  ex- 

Cp  = - j — — —  (9)  cept  near  the  trailing  edge.  The  predicted  results  for  the 

2  1 -stator/  1-rotor  configuration  show  larger  pressure  fluc- 

.  -  tuations  than  for  the  3-stator/4-rotor  configuration  and 

Figure  3  illustrates  the  predicted  and  experimental  the  experimental  data.  The  larger  fluctuations  are  due 
values  of  the  unsteady  pressure  amplitude  coefficient  to  pressure  signals  which  do  not  decay  for  even  blade 
for  the  surface  of  the  stator.  Good  agreement  exists  count  configurations  [4]  and  the  reflective  exit  bound- 
between  the  predicted  3-stator/4-rotor  results  and  the  ary  conditions  used  in  the  numerical  procedure.  Fig- 
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ure  4  illustrates  the  predicted  and  experimental  values 
of  the  unsteady  pressure  amplitude  coefficient  for  the 
rotor.  Good  agreement  exists  between  the  predicted  3- 
stator/4-rotor  results  and  the  experimental  data  except 
near  the  pressure  surface  trailing  edge,  where  the  pres¬ 
sure  variations  are  greater  than  those  observed  experi¬ 
mentally.  Again,  the  predicted  pressure  fluctuations  for 
the  1-stator/ 1-rotor  configuration  are  larger  than  those 
predicted  in  the  3-stator/4-rotor  configuration  as  well 
as  the  experimental  data.  Figures  1-4  illustrate  that 
the  stator/rotor  blade  count  ratio  has  little  effect  on 
the  time-averaged  pressure  field,  but  has  considerable 
impact  on  the  unsteady  pressure  field. 

Figure  5  shows  the  predicted  3-stator/4- rotor  and 
1-stator/ 1-rotor  time-averaged  temperature  coefficient 
distributions  compared  with  the  experimental  data  re¬ 
ported  by  Butler  et  al  [1]  for  the  rotor.  Also  included  in 
Fig.  5  are  the  two-dimensional  numerical  results  of  Rai 
and  Dring  [2]  and  Krouthen  and  Giles  (rescaled)  [5]. 
The  temperature  coefficient,  Ct,  is  defined  as  [2,7]: 


rp _ rr 

Ct  =  r - V  (10) 

where  T  is  the  time-averaged  local  temperature  and 
TavgrU  is  the  time-averaged  mid-span  temperature  at 
the  rotor  leading  edge.  This  definition  of  the  tem¬ 
perature  coefficient  facilitates  comparison  between  the 
present  solutions  and  the  predicted  results  of  Rai  and 
Dring  [2],  Krouthen  and  Giles  [5],  and  the  experimen¬ 
tal  data  [1].  The  time-averaged  temperature  coefficient 
distribution  for  the  1-stator/ 1-rotor  hot  streak  calcula¬ 
tion  is  essentially  the  same  as  that  reported  by  Rai  and 
Dring  [2]  except  for  a  small  difference  on  the  pressure 
surface.  The  discrepancies  on  the  pressure  surface  are 
probably  due  to  the  difference  in  airfoil  scaling  strate¬ 
gies  discussed  eariler.  The  predicted  time-averaged  tem¬ 
perature  coefficient  distribution  for  the  pressure  surface 
of  the  3-stator/4-rotor  simulation  is  similar  to  that  of 
the  1-stator/l-rotor  calculation,  but  the  suction  surface 
temperature  coefficient  distribution  shows  lower  time- 
averaged  temperatures  than  that  predicted  in  the  1- 
stator/ 1-rotor  calculation.  The  lower  suction  surface 
temperatures  predicted  in  the  3-stator/4-rotor  simula¬ 
tion  may  be  the  result  of  introducing  one  hot  streak 
every  third  stator  passage  instead  of  a  hot  streak  in 
each  stator  passage.  The  temperature  coefficient  distri¬ 
bution  obtained  by  Krouthen  and  Giles'[5]  is  similar  to' 
the  present  3-stator/4-rotor  prediction. 

The  segregation  of  the  hot  gases  to  the  pressure  side 
of  the  rotor  passage  has  been,  until  this  time,  believed 


to  be  caused  soley  by  the  relative  inlet  angle  difference 
between  the  hot  streak  and  the  surrounding  fluid  at 
the  inlet  to  the  rotor  passage  [1],  Results  of  the  two- 
dimensional  simulation  confirm  experimental  observa¬ 
tions  which  indicate  that  the  static  pressure,  total  pres¬ 
sure,  and  absolute  flow  angle  are  the  same  in  the  hot 
streak  as  in  the  surrounding  fluid  at  the  stator  exit. 
The  difference  in  the  total  temperature  between  the  hot 
streak  and  the  surrounding  fluid  results  in  an  absolute 
velocity  difference  equal  to  the  square  root  of  the  ratio 
of  the  hot  streak  and  surrounding  fluid  temperatures.  In 
the  relative  frame,  a  difference  in  the  rotor  relative  inlet 
angle  and  velocity  results.  As  a  result,  the  hot  streak 
fluid  moves  towards  the  pressure  surface  at  a  higher 
relative  velocity  compared  to  the  surrounding  fluid  [1]. 
Time-averaged  static  temperature  contours  for  the  rotor 
passage  of  the  two-dimensional  3-stator/4-rotor  config¬ 
uration  are  shown  in  Fig.  6.  Although  the  hot  streak 
fluid  accumulates  near  the  pressure  side  of  the  rotor, 
the  hot  fluid  does  not  penetrate  the  boundary  layer  to 
the  surface  of  the  rotor  blade.  The  thin  layer  of  cooler 
fluid  which  exists  between  the  hot  fluid  and  rotor  pres¬ 
sure  surface  helps  explain  the  flat  temperature  coeffi¬ 
cient  profiles  shown  in  Fig.  5  for  the  two-dimensional 
solutions. 


All  of  the  two-dimensional  numerical  simulations 
shown  in  Fig.  5  predict  nearly  equal  time-averaged  tem¬ 
peratures  on  the  pressure  and  suction  surfaces  and  fail 
to  predict  the  increased  rotor  pressure  surface  tem¬ 
peratures  observed  experimentally.  Clearly,  an  addi¬ 
tional  mechanism  other  than  the  rotor  incidence  vari¬ 
ation  through  the  hot  streak  is  accounting  for  the  ex¬ 
perimentally  observed  temperature  rise  on  the  pressure 
side  of  the  airfoil. 


To  gain  better  insight  into  the  unsteady  phenomena 
within  the  turbine  stage,  animated  sequences  of  the  flow 
field  were  created.  Figure  7  shows  the  temperature  field 
at  one  instant  in  time  during  the  animated  sequence. 
The  hot  streak  is  undisturbed  as  it  migrates  through 
the  stator  passage,  except  that  the  width  of  the  hot 
streak  decreases  due  to  flow  acceleration.  The  geometry 
of  the  turbine  rotor  blades  is  such  that  the  hot  streak 
impinges  first  on  the  suction  surface,  then  wraps  around 
the  leading  edge  and  moves  along  the  pressure  surface. 
As  a  result  of  this  movement,  the  rotor  pressure  surface 
temperature  does  not  peak  until  some  time  after  the 
rotor  blade  interacts  with  the  hot  streak.  The  hot  streak 
assumes  a  V  shape  as  it  continues  to  convect  through 
the  rotor  passage. 
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Figure  5:  Predicted  and  experimental  time-averaged 
surface  temperature  for  rotor  (2-D) 


Figure  6:  Time  averaged  static  temperature  contours 
for  2-D  hot  streak  (3-stator/4-rotor) 


Three-Dimensional  Simulation 

A  1-stator/ 1-rotor  three-dimensional  hot  streak  simu¬ 
lation  with  a  hot  streak  temperature  20  percent  greater 
than  the  free  stream  temperature  has  also  been  per¬ 
formed  to  establish  the  three-dimensional  hot  streak  mi¬ 


Figure  7:  Static  temperature  contours  for  2-D  hot  streak 


gration  effects  on  the  time-averaged  rotor  surface  tem¬ 
perature  distribution.  The  three-dimensional  calcula¬ 
tion  was  performed  on  the  NAS  Cray  2  supercomputer 
located  at  the  NASA  Ames  Research  Center.  This  cal¬ 
culation  required  .000263  seconds  per  grid  point  per  it¬ 
eration  computation  time.  Seven  cycles  at  2000  itera¬ 
tions  per  cycle  were  needed  to  obtain  a  time-periodic 
solution. 

For  the  three-dimensional  simulation,  the  stator  grid 
system  was  constructed  with  101  x  21  grid  points  in 
each  spanwise  O-grid  and  58  x  31  grid  points  in  each 
spanwise  H-grid.  The  rotor  grid  system  was  constructed 
with  101  x  21  grid  points  in  each  spanwise  O-grid  and 
60  x  31  grid  points  in  each  spanwise  H-grid.  A  total 
of  51  O-H  grid  planes  were  distributed  in  the  spanwise 
direction.  The  rotor  region  had  a  tip  clearance  grid 
system  that  contained  101  x  11  grid  points  in  each  of 
7  spanwise  locations.  A  total  of  410,677  grid  points 
were  used  in  the  three-dimensional  simulation.  A  wall 
spacing  of  5.0  x  10-4  inches  was  used  in  the  blade-to- 
blade  direction,  while  a  wall  spacing  of  3  0  x  10~2  inches 
was  used  in  the  spanwise  direction.  The  experimental 
rotor  has  a  tip  gap  equal  to  1%  of  the  rotor  span,  while 
the  tip  clearance  in  the  three-dimensional  simulation  is 
equal  to  3.9%  of  the  rotor  span.  Illustrations  of  the 
three-dimensional  computational  grid  topology  can  be 
found  in  Refs.  [3,15]. 

As  part  of  the  numerical  investigation,  the  predicted 
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Figure  8:  Time-averaged  pressure  coefficient  distribu¬ 
tion  for  stator 


unsteady  pressures  and  temperatures  have  been  time- 
averaged  and  compared  to  the  experimental  data  [14]. 
Figure  8  illustrates  the  predicted  and  experimental 
time-averaged  pressure  coefficient  distributions  on  the 
stator  at  2,  12.5,  25,  50,  75,  87.5  and  98%  span  loca¬ 
tions.  In  general,  there  is  good  agreement  between  the 
predicted  results  and  the  experimental  data.  The  pre¬ 
dicted  results  also  agree  well  with  those  presented  by 
Rai  [3]  and  Madavan  et  al  [15]. 

Predicted  and  experimental  time-averaged  pressure 
coefficient  distributions  for  the  rotor  at  the  2,  12.5, 
25,  50,  75,  87.5  and  98%  span  locations  are  shown  in 
Fig.  9.  Good  agreement  between  the  predicted  results 
and  the  experimental  data  is  observed  from  the  hub  to 
the  mid-span  location  of  the  rotor.  Some  discrepancies 
between  the  predicted  and  the  experimental  pressure 
distributions  are  evident  on  the  suction  surface  of  the 
rotor,  however,  from  the  mid-span  location  out  to  the 
tip.  These  discrepancies  are  probably  caused  by  the 
stronger  secondary  flows  which  arose  in  the  calculation 
due  to  the  relatively  large  tip  clearance. 

Figure  10  shows  the  predicted  and  experimental  val¬ 
ues  of  the  unsteady  pressure  amplitude  coefficient  for 


o  -  EXP  DATA  <t>  =  .78  [13] 

—  1-STATOR/ 1- ROTOR  <t>  -  .78 


Figure  9:  Time-averaged  pressure  coefficient  distribu¬ 
tion  for  rotor 


the  mid-span  section  of  the  stator.  The  predicted  suc¬ 
tion  surface  pressure  fluctuations  near  the  trailing  edge 
are  smaller  than  those  of  the  experimental  data  and  the 
two-dimensional  results  shown  in  Fig.  3.  In  addition, 
the  pressure  amplitude  coefficient  on  the  pressure  sur¬ 
face  does  not  rise  near  the  trailing  edge  compared  to 
the  experimental  data  and  predicted  two-dimensional 
results  (Fig.  3).  Figure  11  illustrates  the  predicted  and 
experimental  values  of  the  unsteady  pressure  amplitude 
coefficient  for  the  mid-span  section  of  the  rotor.  Fair 
agreement  exists  between  the  predicted  results  and  ex¬ 
perimental  data  as  the  predicted  results  show  larger 
pressure  fluctuations  than  those  observed  experimen¬ 
tally  and  predicted  two-dimensionally  (Fig.  4).  The 
stator  and  rotor  pressure  amplitude  results  shown  in 
Figs.  10  and  11  are  similar  to  the  coarse  grid  predic¬ 
tions  of  Rai  [3].  Madavan  et  al  [15]  showed  that  the 
pressure  amplitudes  on  both  the  stator  and  rotor  are 
sensitive  to  the  computational  grid  density.  The  results 
shown  in  Figs.  10  and  11  are  predicted  using  the  same 
fine  computational  grid  density  as  Madavan  et  al  [15]. 
However,  the  wall  spacing  used  in  the  current  calcula¬ 
tion  was  increased  in  order  to  yield  a  unif  rm  grid  near 
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Figure  10:  Pressure  amplitude  coefficient  for  mid-span  Figure  12:  Predicted  and  experimental  time- averaged 
stator  section  surface  temperature  for  rotor  mid-span  section 


Figure  11:  Pressure  amplitude  coefficient  for  mid-span 
rotor  section 


mid-span/mid-gap,  which  was  necessary  to  maintain  a 
circular  hot  streak.  Comparison  between  the  current 
results  and  those  of  Rai  [3]  and  Madavanet  al  [15]  show 
that  the  predicted  pressure  amplitude  is  not  only  sensi¬ 
tive  to  computational  grid  density,  but  also  to  the  span- 
wise  distribution  of  grid  points. 

Figure  12  illustrates  the  predicted  and  experimental 
values  of  the  time-averaged  temperature  coefficient  dis¬ 


tribution  for  the  mid-span  section  of  the  rotor.  Also 
included  in  Fig.  12  are  the  three-dimensional  invis- 
cid  (with  viscous  modelling)  results  of  Takahashi  and 
Ni  [6,7].  The  three-dimensional  simulation  correctly 
predicts  that  the  pressure  surface  of  the  rotor  is  sub¬ 
ject  to  much  higher  time-averaged  temperatures  than 
the  suction  surface.  In  addition,  the  present  three- 
dimensional  calculation  correctly  predicts  the  rapid  fall- 
off  in  the  time-averaged  temperature  on  the  suction  sur¬ 
face.  Comparison  between  the  current  results  and  the 
solution  of  Takahashi  and  Ni  shows  that  both  three- 
dimensional  techniques  predict  similar  trends  with  the 
current  Navier-Stokes  technique  providing  more  de¬ 
tailed  agreement  with  the  experimental  data.  This  com¬ 
parison,  along  with  the  two-dimensional  results,  implies 
that  the  interaction  between  the  secondary  and  wall 
layer  flows  in  the  rotor  passage  and  the  hot  streak  is 
an  important  factor  in  the  time-averaged  rotor  surface 
temperature  distribution.  Furthermore,  this  compari¬ 
son  highlights  the  importance  of  accurately  predicting 
the  secondary  and  wall  layer  flows  through  the  turbine 
stage  in  order  to  correctly  predict  the  rotor  surface  tem¬ 
perature  distributions. 

The  predicted  rotor  surface  temperature  coefficient 
contours  are  illustrated  in  Fig.  13,  while  the  experi¬ 
mental-  contours  [1]  are  shown  in  Fig.  14.  The  pre¬ 
dicted  contour  patterns  are  similar  to  the  experimen¬ 
tal  contours,  except  near  the  pressure  surface  trailing 
edge  where  the  numerical  contours  do  not  close  as  the 
experimental  data  indicates.  This  difference  between 
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Figure  13:  Predicted  time-averaged  surface  temperature  contours  for  rotor  airfoil 
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Figure  14:  Experimental  time-averaged  surface  temperature  contours  for  rotor  airfoil  [1] 


the  prediction  and  the  experiment  coincides  with  that 
shown  in  Fig.  12  for  the  mid-span  pressure  surface  near 
the  trailing  edge.  This  discrepancy  is  probably  due  to 
the  difference  in  flow  coefficients  between  the  simula¬ 
tion  (<£  =  .78)  and  the  experiment  ( <j>  =  .68),  resulting 
in  a  rotor  inlet  flow  with  5  degrees  more  positive  inci¬ 
dence  in  the  simulation  than  in  the  experiment.  Both 
the  numerical  and  experimental  results  illustrate  that 
the  hot  fluid  spreads  over  the  entire  pressure  surface  of 
the  rotor,  while  on  the  suction  surface  the  hot  fluid  is 
generally  confined  to  the  mid-span  region  of  the  rotor 
surface. 

As  was  performed  for  the  two-dimensional  simula¬ 
tion,  animated  sequences  of  the  flow  field  were  created 
from  results  of  the  three-dimensional  simulation.  A  se¬ 
quence  from  the  animation  of  the  T  =1.05  temperature 
isotherm  at  four  instants  in  time  during  a  global  cycle 
are  shown  in  Figs.  15  and  16.  These  four  instants  in 
time  correspond  to  0,  25,  50  and  75%  of  the  global  cy¬ 
cle,  where  the  global  cycle  is  equal  to  the  rotor  blade 
moving  a  distance  equal  to  one  stator  pitch.  Figure  15 
shows  the  isotherm  from  a  viewing  position  that  dis¬ 
plays  the  pressure  surface  of  the  stator  blades  and  the 
suction  surface  of  the  rotor  blades.  Figure  16  shows  the 
isotherm  from  a  viewing  position  that  displays  the  suc¬ 


tion  surface  of  the  stator  blades  and  the  pressure  surface 
of  the  rotor  blades.  This  sequence  of  pictures  illustrates 
the  migration  of  the  hot  streak  through  the  stator  pas¬ 
sage  and  bow  it  is  broken  into  discrete  spherical  eddies 
as  it  interacts  with  the  passing  rotor  blades.  The  hot 
fluid  remains  on  the  pressure  surface  of  the  rotor  for  a 
long  period  of  time  after  it  encoun  ers  the  hot  streak, 
eventually  migrating  towards  the  pressure  surface  tip, 
where  it  leaks  over  onto  the  suction  surface. 

Conclusions 

Experimental  studies  have  shown  that  combustor  hot 
streaks  can  significantly  affect  the  secondary  flows  and 
wall  temperature  of  the  first  stage  turbine  rotor.  To 
gain  insight  into  the  secondary  flow  and  heat  transfer 
effects  on  a  first  stage  rotor  due  to  a  combustor  hot 
streak,  numerical  simulations  of  hot  streak  migration 
through  a  turbine  stage  have  been  performed.  The  goal 
of  these  simulations  has  been  to  understand  the  physical 
mechanisms  which  control  the  migration  of  hot  streaks 
through  a  turbine  stage  and  lead  to  “hot  spots”  on  the 
rotor  pressure  surface. 

Results  have  been  presented  for  both  two-dimensional 
and  three-dimensional  simulations  of  hot  streak  migra¬ 
tion  through  a  turbine  stage.  The  predicted  results  from 


Figure  16:  T  —  1.05  Isotherm  for  3-D  hot  streak  ( rotor  pressure  side) 


the  two-dimensional  simulations  indicate  that  blade 
count  ratio  has  little  effect  on  predicted  time-averaged 
surface  pressure  and  temperature  distributions,  but  a 
substantial  effect  on  the  unsteady  flow  characteristics. 
These  two-dimensional  simulations  fail  to  predict  the 
experimentally  observed  increase  in  time-averaged  pres¬ 
sure  surface  temperatures  due  to  the  hot  streak.  Results 
of  the  three-dimensional  hot  streak  simulation,  however, 
confirm  experimental  observations  that  high  tempera¬ 
ture  hot  streak  fluid  accumulates  on  the  pressure  surface 
of  the  rotor  blades,  resulting  in  a  high  time-averaged 
surface  temperature  “hot  spot”.  The  inability  of  the 
two-dimensional  simulations  to  predict  the  increased  ro¬ 
tor  pressure  surface  temperatures  underscores  the  im¬ 
portance  of  including  three-dimensional  viscous  effects. 
In  addition,  the  influence  of  secondary  flows  on  the  pre¬ 
dicted  three-dimensional  migration  pattern  of  the  hot 
streak  has  been  shown  to  be  consistent  with  experimen¬ 
tal  data. 
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